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Abstract. In generic Hamiltonian systems with a mixed phase space chaotic 
transport may be directed and ballistic rather than diffusive. We investigate one 
particular model showing this behaviour, namely a spatially periodic billiard chain in 
which electrons move under the influence of a perpendicular magnetic field. We analyze 
the phase-space structure and derive an explicit expression for the chaotic transport 
velocity. Unlike previous studies of directed chaos our model has a parameter regime 
in which the dispersion of an ensemble of chaotic trajectories around its moving center 
of mass is essentially diffusive. We explain how in this limit the deterministic chaos 
reduces to a biased random walk in a billiard with a rough surface. The diffusion 
constant for this simplified model is calculated analytically. 



1. Introduction 

With the term directed chaos we refer to extended chaotic systems in which the time- 
averaged velocity of almost all chaotic trajectories approaches for long time a non-zero 
constant v c ^ 7^ 0. In contrast to the well-known deterministic chaotic diffusion, the 
chaotic transport is ballistic and directed in such systems, (x) = v^t. For ID driven 
Hamiltonian systems it has been shown that directed chaos may exist if the periodic 
driving is such that (i) all generalized time-reversal symmetries are broken [11 and (ii) 
the system has a mixed phase space in which regular and chaotic dynamics coexist 

Examples of this type have been investigated in a number of recent publications 
[H El El Ej and closely related phenomena were realized experimentally with cold atoms 
in pulsed optical potentials EH] - 

In the present paper we address directed chaos in a different situation which might 
be relevant in the context of solid-state physics. We study quasi one-dimensional billiard 
chains with (i) a transversal magnetic field breaking time-reversal symmetry and (ii) 
an asymmetric configuration of scatterers which leads to a mixed phase space. We 
will show that in these systems the nature of transport is entirely different for the 
forward and the backward direction. In one of the two directions transport is due to 
regular orbits skipping along a boundary of the billiard chain (waveguide) without back 
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scattering. In the other direction there is strong back scattering and the dynamics is 
chaotic. Consequently all dynamical properties such as the average transport velocity 
or the superimposed spreading of a distribution of particles are different for the two 
transport directions. It is this special property which makes the systems we consider 
interesting and potentially useful for controlling nanoscale electronic transport. 

We can considerably extend a preceeding study of directed chaos in magnetic 
billiards [7\ because the specific geometry which we propose leads to a particularly simple 
phase-space structure. This allows for detailed analytical calculations and provides 
sufficient insight in order to control the transport properties of our model with a few 
geometric parameters. For example, we obtain an explicit expression for the chaotic 
transport velocity. We can also understand in some detail the velocity dispersion of the 
chaotic trajectories. In known examples for directed chaos this dispersion was 

found to be non-Gaussian, with clear signatures of Levy walks. In contrast, we identify 
a parameter that can be tuned such that the dispersion becomes essentially diffusive 
while t> c h 7^ remains constant. In other words, there really is a clear separation of 
scales between directed transport, (x) ~ t, and undirected broadening of an ensemble, 
(Ax) ~ t 1 / 2 . This might be a desirable feature in applications. In this regime of biased 
diffusion the chaotic dynamics can be approximated by a non-deterministic random walk 
in a billiard with a rough boundary, although our original model is deterministic and 
has no disorder. Rough billiards are nowadays standard models for electron dynamics 
on mesoscopic scales [HI El - However, the system which we study here seems to be the 
first such model with a mixed regular-random phase space and is thus interesting also 
in its own right. 

In the following two main sections we introduce and analyze our billiard models 
with directed chaos. Sec. |2 is devoted to a periodic billiard chain with transversal 
magnetic field in which the dynamics is deterministic, while in Sec. El we study a billiard 
with a rough surface and non-deterministic dynamics. The connection between these 
two models is explained in Sec. 12.31 Sec. H] closes our paper with a summary and some 
concluding remarks. 

2. Periodic chain of magnetic Sinai billiards 

In Fig. ^ we show an example for an asymmetric magnetic billiard which illustrates 
qualitatively the effect we investigate. The system consists of a waveguide in which 
one wall is perfectly straight, yi(x) = b. The opposite wall is distorted periodically, 
V = 3/2 (z) > 0, with 



Specifically in our case this is achieved by a chain of semi-circular obstacles with radius 
R = Vxaax.1 but other geometries may lead to similar results^. 

| For example, we found similar behavior when rectangular unit cells are extended by attaching on 
one side a semicircle instead of cutting it out . 




(1) 
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Figure 1. In a waveguide with obstacles placed only along one of the walls a 
transversal magnetic field leads to the coexistence of regular and chaotic trajectories. 
Two corresponding examples are shown with a dashed and a full line, respectively. Note 
that for the same magnetic field and particle charge the trajectories can be traversed 
only with the indicated orientation. The relevant parameters of the model are a = 1. 2b 
(spatial period), R = 0.56 (radius of obstacles) and r — 2b (cyclotron radius). 



Energy is conserved, and we will consider the dynamics on the energy shell 
corresponding to the velocity Vq. For numerical calculations we use dimensionless units 
in which vq = b = 1, i.e., we measure velocity in units of vq, lengths in units of b and 
time in units of b/vQ. As a consequence, the geometric length of a trajectory is equal 
to the elapsed time. Due to a perpendicular magnetic field the trajectories of electrons 
in the waveguide consist of circular arcs with cyclotron radius r = mvo/eB. In the 
following we will use the cyclotron radius to parameterize the magnetic field. 

The magnetic field leads to a special set of trajectories skipping along the clean wall 
of the channel (dashed line in Fig. [TJ. It is a matter of simple geometry to describe the 
transport due to these regular trajectories which by our convention of the magnetic field 
is directed to the left (negative transport velocity). In contrast, trajectories colliding 
with the distorted wall of the waveguide can be chaotic such as the example shown with a 
full line in Fig.[TJ The figure suggests that such trajectories are transporting to the right, 
that is opposite to the skipping ones. We will see below that this is indeed the case. 
This system shows directed chaos, and the transport due to the chaotic trajectories 
compensates for the regular skipping orbits such that no net transport results if one 
averages over all possible initial conditions. 

2.1. Phase space structure 

A phase-space point on the energy shell is completely described by the position (x, y) 
in the waveguide and the angle <p between the velocity vector and the x-axis. In these 
coordinates, the invariant measure is 

d/i = dx dy dip . (2) 

In order to obtain a general impression of the dynamical scenarios supported by our 
model we show a Poincare section (Fig. |2J) obtained by marking y and (p whenever a 
trajectory crosses the boundary of a unit cell, x = mod a, i.e., whenever x is half 
way between two obstacles. The two panels correspond to crossings with positive and 
negative transport velocity 

v x = v cos (p . (3) 
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Figure 2. Phase space portrait of the system of Fig. 2] y and tp are shown for 
x = Omoda and (a) cosf > 0, (b) cos tp < 0. The data represent one chaotic and four 
skipping trajectories. 



Note that the invariant measure of the resulting Poincare map is 

d/ij_ = dy dip cos if . (4) 

In the Poincare section we observe a big chaotic component which contributes to both 
directions. In fact the randomly scattered points in Fig. and b are the trace of a 
single chaotic trajectory, namely the one shown in Fig. ^ In addition we see in Fig. ^Bp 
also the traces of a number of skipping orbits. These trajectories have a second constant 
of the motion besides the energy, namely the y-component of the current arc's center 

Vc(v,<p) = y-rcosip. (5) 

Therefore the skipping orbits constitute a regular (integrable) component of phase space 
which comprises all phase-space points satisfying 

cos (p < 1 , (bj 

r 

that is all arcs whose lowest point is still above the top of the obstacles, y c (y, <f>) — r > 
l/max- It is clear that trajectories with such initial conditions can never be scattered. One 
may argue that there can exist also skipping orbits violating the condition in Eq. ((BJ). 
Say the lowest point of the current arc is below y max but falls into the gap between 
two obstacles. However, if one continues the trajectory it will finally hit one of the 
scatterers, at least if the longitudinal extension As of the arcs is not rationally related 
to the spatial period a. Thus, exceptional skipping orbits within the chaotic component 
may exist but must have measure zero in phase space. 

Are all orbits either chaotic or skipping? There is no reason to expect this. 
Typically, the boundary between a regular and a chaotic component in phase space 
is fractal, with hierarchies of smaller and smaller stable islands embedded into the 
chaotic sea. Also deep within the chaotic sea stable periodic orbits with small islands 
surrounding them can exist. Such islands are clearly visible for relatively small cyclotron 
radius, e. g. r = 1. However, for the parameters used in Figs. ^ and EJ we were unable 
to detect stable islands even after increasing the resolution of the phase space portrait 
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considerably. For growing cyclotron radius it must be expected that stable islands are 
less and less important, as for zero magnetic field (r = oo) the system reduces to the 
Sinai billiard (Lorentz gas) which is known to be fully hyperbolic [H]. Hence, we will 
neglect the influence of additional stable islands and assume that for sufficiently large 
r the phase space of our model has the simple structure which is visible in Fig. 121 All 
analytical calculations below will be based on this assumption. 



2.2. Average chaotic transport velocity 

We will now calculate the asymptotic chaotic transport velocity. For this purpose it 
is useful to restrict the dynamics to a unit cell by taking x modulo a. After this 
transformation the phase-space volume of the system is finite but the transport velocity 
Vx(t) = Vo cos <p(t) is unchanged. Assuming ergodicity within the chaotic component of 
the unit cell, the long-time average over the transport velocity converges for almost all 
chaotic trajectories to the corresponding phase-space average, 

v ch = limr 1 f o\t'v x (t') 

= v (cosip) ch . (7) 

This is illustrated in Fig. El As common in systems with a mixed phase space, the 
assumption of ergodicity within the chaotic component cannot be justified analytically. 
Numerical evidence in favour of it has been gathered for our model by checking that 
individual chaotic trajectories cover the available phase-space region uniformly. 

In general there is no reason to expect v c ^ = in Eq. (JJJ). However, an average over 
cos ip is obviously zero, if for any point (x,y,(p) in the domain of integration also the 
point (x, y, ip + ir) with opposite velocity direction contributes. For example this is the 
case if the integration extends over the entire phase space of the unit cell, 

(cosy?) uc = 0. (8) 

Hence, in a completely chaotic billiard we would have t> cn = 0, i.e., there is no directed 
chaos. Similarly, in any system with a single chaotic component and time-reversal 
symmetry there can be no directed chaos even if stable islands are present. On the other 
hand directed chaos can exist in fully chaotic systems with time-reversal symmetry if 
the phase space contains two invariant chaotic components, see JE] for an example. 

In our case we can decompose the full phase space into the chaotic component and 
the skipping orbits. Thus we have 

fi uc (cos<p) uc = fJ rcg (cosy?) rcg + fi c h(cosy?) ch (9) 

and 

^uc ^reg ~t~ ^ch 5 (10) 

where f2j = /jd/x denotes the phase-space volume of the set i. Substitution into Eq. (j7j) 
yields 

Vt^„ (COS (£>)rc,„ 
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Figure 3. For the parameters of Fig.^the time-averaged velocity of 5 different chaotic 
trajectories is shown as a function of time. For all trajectories the average velocity 
approaches the value predicted by Eq. (|17J) which is indicated by the horizontal line. 



which is the analogue to the sum rule derived in j2] for the chaotic transport velocity in 
driven ID systems. However, in the present case we are able to calculate explicitly. 
First we note that the phase-space volume of the unit cell is 



27T (ab-^R 2 ) . 



(12) 



Next we introduce the characteristic function of the regular component which, according 
to Eq. (JBJ), is given by 

'y — R — r 



Xres(x, y,(f) = © 



cos <p 



(13) 



The volume of the regular component is then 



n 



d/x x re g(£, y, <p) 



reg J AregV 

a i d<p / d y e{ 











r 



COS if 



j dep (b — [R + r + r cos <£>]) 

J n—9 



where 



2ar(sm6 — 9 cosi 
b-R\ 



9 = arccos 1 — 



(14) 



(15) 



denotes the maximum deviation of the limiting skipping orbit from the horizontal. 
Similarly we have 



^reg( COS( ^)i 



dnXvcg(x,y,Lp) COS if 



l-TT + V 

. / dip {b — [R + r + r cos if]) cos ip 

J-k-6 
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= — ar {9 — sin 9 cos 9) . (16) 

Substitution of these expressions into Eq. (|TT|) finally yields 

Vq 9 — sin 9 cos 9 

^ch = t \ (17) 

2 (Tr/ar) ( a & _ _ ( sin g _ g cos 

for the chaotic transport velocity This result is confirmed numerically in Figs. El and |3J 
For small magnetic field, when our assumption of a two-component phase space is 
justified best, we can approximate Eq. (|T7jl by the leading order in the cyclotron radius. 
We find 



e - (is) 

and 

4a l2(b-R) 3 , . 

This shows explicitly that the chaotic transport vanishes in the absence of a magnetic 
field (r — > oo). 

5. Dispersion of chaotic trajectories 

With the above results we can precisely predict the average transport velocity of many 
chaotic trajectories at a given moment in time, or the long-time average over a typical 
trajectory. However, we have no precise information about the dispersion of an ensemble, 
or the rate of convergence of the transport velocity to its asymptotic value. In systems 
with sufficiently fast decay of correlations, the dispersion of chaotic trajectories is 
diffusive, 

([x(t)-v ch t] 2 ) = Dt. (20) 
In this case the distribution of time-averaged velocities 

vt = X(t) = X{0) (21) 
is for long times a Gaussian with variance D/t, 



P(vt) = \/^75 ex P (-[*/2I>] - ^ch] 2 ) • (22) 

In Fig. 13] we compare this hypothesis to the chaotic dispersion for the parameters 
of Fig. n Clear deviations from diffusive behavior are observed in the tails of the 
distribution and for its higher moments. After iterating iV = 10 6 trajectories to t — 10 3 
the average chaotic velocity agrees very well with the value predicted by Eq. (|17|) . 
However, the distribution of time-averaged velocities is asymmetric around t> c h and, 
compared to a Gaussian, large deviations from the mean value are enhanced. For the 
considered time interval the variance of the distribution (+) is fitted well by the diffusive 
prediction D/t. However, with increasing n the higher moments ((v — f C h) n ) of the 
velocity distribution deviate more and more from the corresponding time-dependent 
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Figure 4. For the model of Fig. ^ the dispersion of chaotic transport velocities 
is analyzed, (a) N = 10 6 trajectories have been iterated up to t = 10 3 . The 
resulting distribution of average velocities is displayed as a histogram and compared to 
a Gaussian with the same mean and variance (dotted line) . The prediction of Eq. I|17|) 
for the average velocity is indicated by a vertical line. The numerical mean value 
of the distribution differs from that by 1.5 x 10~ 5 . (b) shows some higher moments 
((TJ — v c h) n ) of the distribution as a function of time. The straight lines indicate the 
corresponding moments of a diffusively spreading Gaussian with D = 1.2. 
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Figure 5. Same as Fig. 0] for R — 0.016 and a = 2R (neighbouring semicircles are 
touching). The fitted diffusion constant is D = 1.2 



Gaussian (x, * vs full lines). This is no surprise as in systems with a mixed phase 
space anomalous diffusion is typical. Indeed, in previous studies of directed chaos in 
such systems the observed velocity distributions were always non-Gaussian [U |21 El E] ■ 
They displayed traces of Levy walks and anomalous diffusion which were even stronger 
than those in Fig. EJ and it was not possible to adjust the model parameters such that 
directed transport coexists with normal diffusion. 

Therefore it is quite remarkable that in our model a slightly modified geometry 
leads to a velocity distribution which numerically could not be distinguished from a 
Gaussian. In Fig. El the analysis of Fig. 0] is repeated for a system in which (i) the 
semicircular obstacles are much smaller than in Fig. d and (ii) neighbouring scatterers 
are tangent rather than separated by a gap (see Fig. [7^,). 

This latter modification removes marginally stable periodic orbits which hit the 
lower wall within the straight sections of the boundary only (Fig. IH}. Up to a certain 
limit it is possible to translate such orbits along the channel without changing their 
shape. Nevertheless they form a set of measure zero in phase space since the angle of 
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Figure 6. The full lines show two marginally stable periodic orbits (one of them is 
periodic up to translation by one unit-cell length). The dotted line is a non-periodic 
chaotic trajectory which remains in the vicinity of a marginally stable orbit for a very 
long period of time. 

the trajectory cannot be varied without destroying the orbit (dotted line in Fig. EJ). It 
is known that marginally stable trajectories are a source of non-exponential decay of 
correlations. They lead to anomalous diffusion even in completely chaotic systems like 
the Sinai or the stadium billiard. In our system generic chaotic trajectories can remain 
in their vicinity for a long time and during these episodes keep the transport velocity 
substantially above or below the average (e. g. v rs 0.5 for the dotted line in Fig. Awhile 
t> c h = 0.12). Therefore it is quite natural that the chaotic dispersion shows much less 
anomalous diffusion when the specific geometry prevents the existence of marginally 
stable trajectories. 

In contrast, it is not immediately clear why reducing the radius of the scatterers 
should further suppress anomalous diffusion in our model. In order to find an intuitive 
explanation for this numerical fact we consider now a different Poincare section defined 
by the intersections of a trajectory with the line y = R. Note that regular skipping 
orbits do not at all intersect this line. In analogy to Eq. (J3J the invariant measure 
restricted to the surface of section is now 

d/i|| = dxd(p simp = dxdc (23) 

with c = cos if. Hence, a typical chaotic trajectory is represented in this Poincare section 
by points (£„,c„) which uniformly cover the area [0,1] x [—1,1]. £ n = (x n /a)modl 
denotes here the position of the intersection relative to the unit cell. Following a 
transition across y = R from below, (£ n _i,c n _i), there will always be an intersection 
from above which is denoted by (£' n , c' n ). Geometry requires c' n = c n _i. 

For small obstacle radius transport is mainly due to the long segments of the 
trajectory above y = R (Fig. Uh)- Each segment is either a single circular arc 
or a combination of two arcs which are identical up to a reflection. Segment n is 
entirely characterized by the value of c n . Therefore the statistics of the c n determine 
the transport properties. In particular, the velocity dispersion will be diffusive if 
correlations between consecutive angles decay sufficiently fast. This argument will be 
made more explicit in Sec. El below. Fig. [7^ shows that for sufficiently small obstacle 
radius (c c n ) indeed decays very fast as a function of n (exponentially fast to a good 
approximation). This observation can be understood as follows: As R — > 0, a segment 
(£ n _i,c n _i) — > (C n > c 'n) traverses more and more unit cells between its end points. As a 
consequence, one expects an extremely rapid variation of £' n as a function of c n _i. This 
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Figure 7. (a) A piece of a chaotic trajectory is shown for R = 0.1b and a = 2R. 
(b) One particular scattering of this trajectory from the lower billiard boundary is 
magnified. Note that on this scale the curvature of the trajectory due to the magnetic 
field is hardly visible, (c) The normalized correlation between consecutive intersection 
angles of the segments with the line y = R is shown with stars for the model of (a) , with 
empty squares for R — 0.016 and and with black dots for the same R but randomized 
(see text). The straight line shows the fitted exponential exp(— 1.65n). 

suggests to replace £' n by a random variable. Fig.[?t: confirms that this approximation is 
very good for the parameters of Fig. or even smaller obstacle radius R. Note that the 
argument leading to this randomization does not depend on the value of the cyclotron 
radius r. In particular, it remains valid also in the absence of a magnetic field, r = oo. 

Next we consider the dynamics in the vicinity of the lower boundary, i.e., for y < R. 
For small obstacle radius it is practically identical to that for r = oo since the ratio 2R/r 
between the size of the unit cell and the (constant) radius of curvature of the trajectory 
vanishes as R — > (Fig. 03). Hence the resulting mapping (£„, c„) — > (£ f n , cf n ) may be 
approximated by a Sinai billiard without magnetic field. The latter is known to be 
hyperbolic jTT] and shows exponential mixing if no marginally stable orbits are present. 
As the correlations should not be enhanced by the additional randomization of £ we 
finally understand why (coc n ) decays exponentially for small obstacle radius R <C r, 
leading to a diffusive velocity dispersion. 

3. A waveguide with one rough boundary 

The results of the previous section suggest that for small obstacle radius R the detailed 
dynamics inside the lower boundary layer of the waveguide is irrelevant. This holds 
at least as long as the dynamics therein is chaotic. In this case the essential physical 
properties of our model are not affected if we replace the boundary layer by an idealized 
rough surface and describe the scattering at the lower wall probabilistically. In the 
present section we follow this approach and obtain in this way a better analytical 
understanding. For the sake of simplicity we will continue to refer to trajectories 
which hit the lower wall as chaotic although the dynamics is not deterministic anymore. 
This is justified since we have just explained how random scattering may arise from 
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sin(cp) sin(<p) 

Figure 8. (a) Typical electron trajectories (dashed line: regular, full line: chaotic) 
for a waveguide with a perpendicular magnetic field (r = 55/3). The lower wall is 
perfectly rough (a = 0) such that the angles before and after a reflection are completely 
uncorrelated. (b) and (c) show a phase space portrait of this system for r = b. 



a deterministic chaotic system in a certain limit. Nevertheless there are also other 
mechanisms such as thermal or quantum fluctuations which may lead to the same 
probabilistic model. 

In order to define our simplified model we must prescribe the probability for a 
trajectory to leave the lower boundary at an angle < <fi < it. Consistency with the 
invariant measure rf/xy, Eq. (J23|) . requires 

Pffl = ^ (0 < < tt) . (24) 

In terms of c = cos <fi this amounts to 

P(c) = l (-I<c<+1) (25) 

for the probability density averaged over all chaotic trajectories. We account for possible 
correlations between the segments of a trajectory by introducing a probability a for a 
specular reflection at the lower boundary. The probability density for c n is then 

1 - a 

P(c n ) = a S(c n - c n _!) + . (26) 

While for a = there are no correlations between consecutive values of c n 

r+l 



^ dcP(c)c 



... (a = 0) , (27) 

the possibility of specular reflections leads to exponentially decaying correlations 

(c c„) = — (a ^ 0) . (28) 
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For each segment of a chaotic trajectory we can express the horizontal distance between 
its endpoints as 

Ax(4>) = 2r(sin - sin 0') . (29) 

On the other hand the total length of a segment is 

AZ(0) = 2r(0-0'). (30) 

Here and 0' denote the angles of intersection with the lower and the upper wall of the 
channel, respectively (see Fig. [S^). In terms of the former, the latter is given by 

arccos(cos0 + -) (cos0 < 1 — -) 

r , (31) 
(cos0>l-£) 

where the second line extends the definition of 0' to those segments of a trajectory which 
do not reach the upper wall. 

With the help of the angle we can give an alternative parameterization of the 
chaotic component of phase space. Instead of (x, y, <p) we can refer to one of its points by 
the coordinates (x, 0, I). As mentioned above, is the inclination angle of the segment 
at its initial point, and < / < A/(0) is the path length measured from there. Formally, 
the transformation between (0, 1) and ((p, y) is given by 

f0-I l<^ 

<P(M=\ " Alu) (32) 

[ -</?(0,A/-/) l> 

|/(0,/) = T cos v?(0, 1) — T COS , (33) 

where the second line of Eq. f!32j) expresses the reflection symmetry of the trajectory 
segments. In terms of and I we can express arbitrary averages over the chaotic phase 
space component as 

Jo 7r d0P(0)/ o A;W d/(...) 

/ c Td0P(0)/ o AW d/ 

f"K rA.l((f)) 

(A/)! 1 / d0P(0) / dl (...), (34) 



(• • .)ch 



10 JO 

where (. . .)^ is the average over the probability density P(0). Eq. will be applied 
in the following sections. 

It is useful to note that Eqs. (J2l?j) and (|5U|) allow to understand the chaotic transport 
as random walk along the x-axis without any reference to the transversal motion. The 
discrete steps of this random walk are Ax(0) and the corresponding time increments 
are A/(0)/vo- Step n is chosen according to the probability density given in Eq. ([26)1 . 

3.1. Average chaotic transport velocity 

We will now give an alternative derivation for the chaotic transport velocity in the case 
R = 0. We can apply Eq. (jMj) to Eq. ((7j) and simplify the result using the identity 

r Al(<fi) 

dl cosip((j),l) = Ax(0) (35) 
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which follows from Eq. (|32j) and should also be obvious from the geometrical meaning 
of Ax and A/. In this way we obtain 



Al(4>) 



v ch = v (Al)^ j Q d0P(0) 

= v D (Al)t 1 £d</>P(<t>)Ax{<l>) 
(Ax)^ 



dl cos<p 



v 



(Al), 



(36) 



This representation is very intuitive from the random-walk point of view: a long 
trajectory consists of iV — > oo steps distributed according to P(0). Therfore the 
distance between the end points of the trajectory is AX = N(Ax) ( j ) while the total time 
increment is AT = N(Al) ( f ) /v . The average longitudinal component of the velocity is 
then AX/ AT which is equivalent to Eq. (}3l)j) . Explicit averaging over yields 



(Ax), 



sin i 



d(f) sin (sin < 
J +1 dc Vl - c 2 - r J 1 ^ dc yjl - [c + b/r] 2 ) 



n r 
~2 ~ 4 







b 




7TH 






(-•) 






r 





r 
2< 



sin 9 cos ( 



(37) 



In agreement with Eq. (fTH|) we have set cos 8 = 1 — b/r in order to obtain the last line. 
For the average length of a segment we have 



(A/) = r / d4> sin 0(0-0') 



+i 



dc arccos(c) — r 



-i 



l-b/r 
-1 



r 

7rr 

2 



7T- 

r 



dc arccos(c + b/r 
b 



7i — 2 arcsin 1 

V r 



A 



b - 2- b - 
r V r , 



(sin^-^cos^) 



(3f 



Combining Eqs. (j36|l - (j38|) we find the expected result, namely Eq. (fPTj) with R = (and 
arbitrary a) 

— sin 9 cos ^ 



V 



(39) 
(40) 



2 (nb/r) - (sm6 - 6 cos6) 

2v [2b . , 
3n V r 

Fig. E3 shows Eq. (j3H|l as a function of the cyclotron radius. We observe that for 
small r the properties of the system change qualitatively. In particular the chaotic 
transport velocity is constant, v c ^ = v /2 for r < 6/2. In this regime a third phase 



Directed chaos in magnetic billiards 



14 




Figure 9. Dependence of the chaotic transport velocity and the diffusion constant 
on the cyclotron radius r for random scattering from the lower wall (a = 0). 



space component of pinned cyclotron orbits appears in addition to regular skipping and 
chaotic trajectories. Hence, the phase-space sum rule for t> c h must be modified and 
Eq. (JT7j) is not valid. In contrast, Eq. (|36|) still applies, but now we have 0' = for all 
<p since no chaotic trajectory reaches the upper wall. Therefore we have 



(Ax)^ = 2r(sin 
1 
2 



1 P+ 1 
2r- I dcos0 sin< 



and finally 



= 2 r > ( 41 ) 
= ixr (42) 



^ch = y • (43) 



3.2. Diffusion constant 

Using the method of the previous subsection we can also calculate the diffusion constant 
D of our simplified model. Starting point is the well-known Green-Kubo relation 
expressing D as integral over the velocity autocorrelation function, 

/+oo 
dt(v x (t)v x (0)). (44) 
-co 

Here v x (t) = v cos ip(t) — t> ch denotes the fluctuation of the transport velocity around its 
mean value and the average (. . .) is taken over all possible trajectories. For the moment 
we will assume a = 0, i.e., no correlations between different segments of a trajectory. In 
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this case it suffices to restrict the integral over t to the segment containing the phase- 
space point (0 O , l ) at t — 0, that is to the interval —l < v t < A/(0 O ) — l . We use 
Eq. (jHIJ) to average over (0 O , ^o), make the substitution I = l + t> t and apply Eq. (fHKJ) . 
In this way we find 

A>=o = / #oP(0o) / d/od/^ 1 (45) 



<AZ) Jo 

X («„ COSV2(0O, - ^ch) («o cos(0 o , lo) - v ch ) 



f d0 o P(0 o ) (Ax(0o)-— A/( 



(A/)^ Jo \ v 

([Ax - (v ch /v )Al} 2 ) (f> 

~ Vo W* ' m 

This result reduces the calculation of the diffusion constant to a simple ID quadrature 
which is numerically straightforward. An explicit result is available for strong magnetic 
field 

^-^ K). <«> 

which reproduces the initially linear dependence of the diffusion constant on the 
cyclotron radius (Fig. |UJ). 

If a 0, several consecutive segments of a trajectory can be identical. Accordingly, 
the section of the trajectory which is correlated with the point (<f) , l ) is longer and 
the integration over I in Eq. (}4"5|) extends over the range — Iq — //A/(0o) < I < 
A/(0 O ) — + vAl((j)o) where //, v > count the identical segments before and after 
the current one. Hence the integral increases by a factor \i + v + 1. For z/ identical 
segments following the current one we need i/ specular reflections from the lower wall 
and finally one random reflection. The probability for this is ex v (1 — a). Similarly, the 
probability for fi identical segments preceeding the current one is a M (l — a). Averaging 
with these probabilities over /i, v yields 

oo 

= \^D . (48) 
1 — a 

In Fig. ^1 we verify these results on the chaotic diffusion constant numerically. 



4. Summary and discussion 

In this paper we have introduced and studied a particular billiard model which shows 
directed chaos. The most appealing feature of this model is the fact that the phase space 
has, to a good approximation, a very simple structure which can be understood both 
intuitively and analytically. It consists of two invariant manifolds, one chaotic and one 
regular. The latter is formed by trajectories skipping periodically along one perfectly 
straight wall of a channel. The other wall has semicircular obstacles or is disordered 
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Figure 10. Some higher moments of the distribution of N = 10 4 time-averaged 
chaotic transport velocities in the magnetic billiard with rough lower wall and r = 2b. 
The probability of specular reflection is (a) a = and (b) a = 3/4. Straight lines 
indicate diffusive spreading with the diffusion constant given by Eqs. I|46|) and JJSJ. 
The latter correctly predicts that D c h is greater by a factor 7 in (b). 



or rough. This leads to strong back scattering and thus to a chaotic or random phase- 
space component. Both, the particular phase-space structure and the fact that time- 
reversal symmetry is broken, are dependent on the presence of a transversal magnetic 
field acting on charged particles such as electrons. The magnetic field should have a 
suitable intermediate strength. For very small magnetic field the phase-space volume 
of the regular skipping trajectories goes to zero and, accordingly, the chaotic transport 
vanishes. For strong magnetic field the phase-space structure changes as more and more 
regular regions appear. Finally, for very strong magnetic field, pinned cyclotron orbits 
cover most of the phase space. In this regime the chaotic trajectories are skipping along 
the lower wall and t> c h approaches a non-zero constant, but at the same time the chaotic 
fraction of phase space goes to zero. 

We have stressed in the introduction the analogies between directed chaos in driven 
systems and in magnetic billiards. However, there exists also an important difference. 
It has been shown that driven Hamiltonian systems with directed chaos can generate 
directed transport in ensembles of particles which have initially a thermal distribution 
in phase space [U El El IH El El • They do so without an external bias force and under 
the influence of a potential which is periodic both in space and in time. This is very 
similar to the concept of stochastic ratchets (Brownian motors) [TT? 1 IT^ l ITK 1 llfij. Hence 
the term Hamiltonian ratchet is appropriate for the models investigated in PHEIEHUE] 
and also for the atom-optics setup studied experimentally in However, magnetic 
billiards with directed chaos are no Hamiltonian ratchets and cannot generate directed 
transport from thermal ensembles. The reason is the conservation of energy. We have 
seen that the total transport from all regular and chaotic phase-space components of 
an energy shell must vanish. Therefore directed transport in magnetic billiards requires 
control over the initial conditions beyond prescribing a certain distribution of energy 
values. For example, electrons could be placed selectively into the chaotic phase-space 
component of our billiard chain if they enter the system from a lead which is attached 
to the lower wall of the waveguide. 
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An interesting extension of the present work would be the inclusion of quantum 
effects. Definitely they play an important role when applications to electronic transport 
in semiconductor nanostructures are the goal. For a given geometry and a given time 
there is always a semiclassical regime in which the effective value of Planck's constant 
is sufficiently small such that the classical results obtained above remain valid. On 
the other hand, for a given finite h there should be interesting deviations from our 
predictions for large times when tunneling between the regular and the chaotic regions 
of phase space may become important and when the presence or absence of disorder in 
the system is crucial. 
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